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Abstract Electromagnetic form factors have long been used to probe the underlying charge and 
magnetization densities of hadrons and nuclei. Traditional three-dimensional Fourier transform 
methods are not rigorously applicable for systems with constituents that move relativistically. 
The use of the transverse charge density is a new, rigorously defined way to analyze electro- 
magnetic form factors of hadrons. This review is concerned with the following issues: what is 
a transverse charge density; how is one extracted one from elastic scattering data; the existing 
results; what is the relationship with other observable quantities; and, future prospects. 
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1 INTRODUCTION 

A truly impressive level of experimental technique, effort and ingenuity has been 
applied to measuring the electromagnetic form factors of the proton, neutron 
(nucleon) and pion (|TJ[2j[3lH]) ■ These quantities are probability amplitudes that 
a given hadron can absorb a specific amount of momentum and remain in the 
ground state, and therefore should supply information about charge and mag- 
netization densities. The text-book interpretation of these form factors is that 
their Fourier transforms are measurements of the charge and magnetization den- 
sities. This interpretation is deeply buried in the thinking of nuclear and particle 
physicists and continues to guide intuition, as it has since the days of the Nobel 
prize-winning work of Hofstadter([2]). Nevertheless, the relativistic motion of the 
constituents of the system causes the text-book interpretation to be incorrect ([6]). 
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The preceding statement leads to several questions, the first being: Is the 

statement correct? If correct, how fast do the constituents actually have to move 
to violate the non-relativistic interpretation? Why is it that the motion of the 
constituents and not that of the entire system that is relevant? The answers to 
these questions are probably displayed within the existing literature. However, 
obtaining general clear answers has been sufficiently difficult that posing even the 
first question of this paragraph would not lead to a unanimous answer. Moreover, 
there is another more important question. If the non-relativistic approach is not 
correct: What is the correct procedure to determine model- independent infor- 
mation regarding hadronic charge and magnetization densities? How should we 
interpret the beautifully precise electromagnetic form factor data being produced 
at Jefferson Laboratory, Bates, Bonn and Mainz? 

The aim of this review is to present and clarify answers to these questions. The 
most immediate question is addressed here. The only way to determine model- 
independent information about the charge distribution is to study the transverse 
charge density p{b) (|7J) which gives the density of a hadron moving at infinite 
momentum at a transverse (to the direction of rapid motion) distance b from the 
transverse center of momentum. It is only in the infinite momentum frame that 
one can define such a center. Information about the charge density in terms of 
longitudinal coordinates is not yet available in a meaningful way. 

1.1 Electromagnetic form factors are not three-dimensional Fourier 
transforms of charge densities 

Physicists are trained to believe that physics is independent of the inertial ref- 
erence frame. In non-relativistic quantum mechanics this invariance allows us 
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to express the wave function as a product of a plane-wave factor (that describes 

the motion of the center of mass) with a function that depends only on internal 

relative coordinates. Thus for two spinless particles: 

*(r l! f 2 ) = e iA ^(r), (1) 

where P is the total momentum, R is the position of the com, and r = r± — r<i- 
Then the non-relativistic form factor Fnr (\q\), which is the probability amplitude 
for the system to absorb a momentum q and remain in its initial state, is given 
by the integral 

FnrM) = I d\\m\ 2 ^ 12 , (2) 



if the masses of the two constituents are the same. One sees that the square of a 
wave function or probability density is being probed. 

However the separation appearing in Eq. (pQ) is not generally valid. For a 
similar but relativistic system the Minkowski-space, Bethe-Salpeter wave func- 
tion &(k,P) (where k is the relative four-momentum and P is the total four- 
momentum) can be written in a compact form if the interaction kernel is given 
by a set of Feynman graphs (8,9), using the Nakanishi integral representation: 

Ht,p) = -4^ r dz r v — g(7,z) r*, (3) 



\/4<k J-l Jo 



7 + m 2 — \M 2 — k 2 — P ■ k z — is 



where m is the mass of each constituent scalar particle and M is the hadronic 
mass. The weight function 3(7, z) itself is not singular, but the singularities of the 
BS amplitude are reproduced by using Eq. ([3]) (jT0[) . The key feature of Eq. ([3]) is 
that the covariant wave function depends explicitly on the total four-momentum 
P. Another way of presenting the contents of Eq. ([3]) involves the relativistic 
boost operator used to obtain the wave function in a moving frame from one in 
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which the hadron is at rest. The Nakanishi representation allows one to perform 

the boost merely by changing the value of P. 

The explicit dependence on P dramatically influences our understanding of 
form factors because the initial and final hadrons have different momentum and 
therefore different wave functions. The presence of different wave functions of the 
initial and final nucleons invalidates a naive probability or density interpretation. 

Assuming that only one of the scalar constituents carry charge, the electromag- 
netic form factor is obtained by evaluating the impulse approximation triangle 
diagram, with the result (fTtJ|) . 

{P + P'YF(Q 2 ) = i J ^(P + P' + kr(k 2 -m 2 )^P-k,P)^P'-k,P'), (4) 
where P' = P + q. 

The form factor of Eq. Q seems to differ markedly from the usual non- 
relativistic expression Eq. ([2]) , expressed in terms of momentum space wave func- 
tions 4>: 

/ii 3 k -, 
^^{k + q/2)m- (5) 

It is comforting that in the non-relativistic limit, defined by the replacement 
1.0 i t.3 1 

+ — , M-2m<2m, (6) 



P° + P 3 2 2m' 
the IMF, light-front wave function becomes identical to the rest frame instant- 
form wave function (jlip . This means © that if the condition Eq. ([6]) is true for a 
given model, then Eq. ([5]) (and therefore Eq. ([2])) is obtained from Eq. @. In this 
case, one may extract p(r) = |</>(r)| 2 from a three-dimensional Fourier transform 
of Fnr. However, the conditions Eq. © are not expected to be obtained for 
hadrons, although they are valid for nuclei ([6|). Furthermore, if the charged par- 
ticles of a given system have different masses, one obtains (in the non-relativistic 
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limit for all constituents) 



F NR (q) = J d 3 r^ Pi (f) e -^/ M , 



(7) 



where Pi,rrii are the charge density and mass of the i'th constituent, and one can 
not obtain the charge density from the form factor. 

1.2 What is a transverse charge density? 

A proper determination of a charge density requires the measurement of a matrix 
element of a density operator. We show here that measurements of the hadronic 
form factor directly involve the three-dimensional charge density of partons in the 
infinite momentum frame, f>oo{ x ~ , b). Before discussing this quantity we need to 
provide a brief introduction to light-cone coordinates. 

In the infinite momentum frame IMF the time coordinate ct = x°/V2 is ex- 
pressed in a frame moving along the negative z direction with a velocity nearly 
that of light using the Lorentz transformation as the variable x + = (x° + x 3 ) / \/2, 
with the usual 7 factor is absorbed by a Jacobean of an integral over volume 
(|12j) . The x + variable is canonically conjugate to the minus-component of the 
momentum operator p~ = (p° — p 3 )/^/2. The longitudinal spatial variable is 
x~ = (x° — x 3 )/\/2 and its conjugate momentum is p + = (p° + p 3 )/\/2. It is this 
plus-component of momentum that is associated with the usual Bjkoren variable. 
The transverse coordinates x, y are written as b with the conjugate p. In the 
literature, these transverse coordinates are sometimes written with a subscript 
_L, but here we shall simply use boldface to denote the transverse components 
of position and momentum vectors. We shall also always use a light-front time 
quantization which sets x + and the plus-component of all spatial variables to zero. 
This means that x~ can be thought of as the longitudinal variable — ^/2x 3 . One 
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extremely useful aspect of using these variables is that Lorentz transformations 

to frames moving with different transverse velocities do not depend on interac- 
tions. These transformations form the kinematic subgroup of the Poincare group 
so that boosts in the transverse direction are accomplished using the transverse 
component of the position operator (as in the non-relativistic theory). 
In the IMF, the electromagnetic charge density J° becomes J + and 

p^x' ,h) = J + (x~,b) = ^e q q(x~ ,b)j + q(x~ ,b) = ^ e q V2q\_(x~, b)q + (x~~, b), 

<? <? 

where q + {x li ) = 7°7 + /v / 2<2 , (^ M ), the independent part of the quark-field operator 
q(x^). The time variable, x + is set to zero. Note the appearance of the absolute 
square of quark field-operators, which is the signature of a true density. 

We are concerned with the relationship between charge density and the elec- 
tromagnetic form factor F(Q 2 ), determined from the current density via 

F(o2) ^pyw.ri (9) 

where the normalization is (p' + , p'|p + , p, ) = 2p + (27r) 3 <5(|/ + — p + )6^(p' — p). 
Spin or helicity dependence is ignored in the present sub-section. We take the 
momentum transfer q a = p' a — p a to be space-like, with the square of the space- 
like four-momentum transfer q 2 = —Q 2 and use the Drell-Yan (DY) frame with 
(q + = 0, Q 2 = q 2 ). No longitudinal momentum is transferred, so that initial and 
final states are related only by kinematic transformations. Moreover, with this 
condition the current operator links Fock-state components with the same number 
of constituents. The matrix element appearing in Eq. © involves the combination 
of creation and destruction operators: b^b — d)d for each flavor of quark, so 
that the valence charge density is probed. The form factor F is independent of 
renormalization scale because the vector current q^q is conserved f)13|) . 
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The spatial structure of a hadron can be examined using states that are trans- 
versely localized ([71fT^ l [T3 |) through a linear superposition 

\ p +,K = O)^AfJ-0^\p+, P ), (10) 

where M is a normalization factor satisfying |A/"| 2 / Jpjfk = Wave packet repre- 
sentations can be used to avoid states normalized to 5 functions (115 pi6p . but this 
leads to the same results as using Eq. (fTOl) . Considering 2p + p~ — p 2 = m 2 > 0, 
(with p~ > 0) one finds p + must be very large because the range of integration 
over p is large. Using such an ultra- large or infinite value of p + , maintains the 
interpretation of a pion moving in the IMF with well defined longitudinal momen- 
tum (|15p . It is in just such an infinite momentum that the parton interpretation 
of a hadron is valid. Setting the transverse center of momentum to zero, Eq. (|10p . 
allows a meaningful transverse distance b. A Fock-space parton representation of 
the position of the transverse center of momentum, R, provides a useful relation 
between longitudinal momentum fractions and transverse positions (|14p : 

R = ^x 4 b 4 , (11) 

i 

where the sum is over the finite number of constituents in a given component. 
Next we relate the charge density 

/ - ^ (p + ,R = 0|poo(a:~,b) |p+,R = 0) 
P °° {X ' b) = (p+,R = 0|p+,R = 0) ' (12) 

to F(Q 2 ). In the DY frame no momentum is transferred in the plus-direction, 
so that information regarding the x~ dependence of the density is not accessi- 
ble. Therefore we integrate over x~ , using the relationship q+(x~ , h)q + (x~ , b) = 
e v + *- e -*P-b 9 t_(0) g+ (0)e*P ,b e- , ^ +:e ~, to find 

p(b) = y'dx-p oo (x-,b) = (p + ,R = 0|p oo (0,b)|p+,R = 0)/(2p + ). (13) 



Transverse Charge Densities 



9 



The use of Eq. Q and the expansion Eq. (|10p leads to the simplification of the 
right-hand-side of the above equation: 



where p(b) is termed the transverse charge density, giving the charge density at 
a transverse position b, irrespective of the value of the longitudinal position or 
momentum. The use of a three-dimensional coordinate space density poo(x~ ,h) 
to obtain the transverse density appeared in Ref. (fTTj) . Previous derivations (UJ 
[T4" l ll3 1[T8|) used a density operator involving transverse position and longitudinal 
momentum variables; see sect. 2. 

If the non-relativistic (NR) limit Eq. Q is valid, one finds that ([6]) 



where p(r) is the square of the wave function obtained from Eq. ([2]). This, as well 
as the integral over the longitudinal coordinate that appears in Eq. illus- 
trates the nature of the transverse density as a reduction of the three-dimensional 
density. Moreover, these integrals appear in the Glauber theory (|19p of high en- 
ergy scattering in which scattering amplitudes are expressed in terms of transverse 
densities. A recent example is found in (20). 

1.3 Simple example 

It is worthwhile to use a simple example (jB) in which a scalar particle is modeled 
as a bound state of two different scalar particles £, x, with a point-like coupling 
such that the interaction Lagrangian is g^f^X- Then the Bethe-Salpeter wave 
function is obtained from Eq. ([3]) by replacing 5(7, z) / \f^n by the simple coupling- 
constant g. The explicit Bethe-Salpeter wave function can then be obtained by 




(14) 




(15) 
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straightforward integration and found to be proportional to the product of Klein- 

Gordan propagators. For this model, the evaluation of the form factor Eq. @ 

can be performed in three different ways. One can simply evaluate the Feynman 

integral, one can take work in the IMF with P 3 — > oo (21) or one can proceed by 

first integrating over k~ . The resulting form factor is the same in all three cases 

([6]). The form factor can be expressed in terms of a three-dimensional integral 

which involves a light-front wave function: ip(x, k) 

1>{X, K) EE 5 [M 2 - ^±^1 - ^±M]"1. (16) 

x 1 — X 

In that case 

l f dx 

F{Q2)= 2(2^ J * K 1 x(1 _ x) V>*(^+(l-x)q)V^,K). (17) 

The relative variables are x, k with x = k + /P + and k = m — «2- The function 
ip(x,n) can be obtained from <I> of Eq. © by integrating over k~ ( 22 , 23| [TT]) . 

One can obtain a relation between the form factor and a coordinate space den- 
sity, by expressing the light-front wave function vp(x, k) in terms of the canonically 
conjugate spatial variable B = b — b2, where b is the transverse position operator 
of the charged constituent. Then ip(x, B) = J e lK ' B ip(x, n), and 

^> = s/^/. 1 ^i** b >i v ^ k (18) 

We express the form factor in terms of the position of the charged constituent 
by using the condition that the center of transverse momentum is set to zero 
(Eq. (jlip with two components): R = = xb + (1 — x)b2 so that B = b/(l — x). 
Expressing the form factor in terms of b gives 

F{Q 2 ) = ]- ( d 2 b ^ . dX |fe — ^— )| 2 e- ib ^. (19) 
Air J Jo x(l-x) 3 l-x n 
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Comparing with Eq. (|14p allows us to identify the transverse density as 



The mass M must be less than 2m for the hadron \& to be stable. We define 
a positive binding fraction e, so that M = 2m — em (0 < e < 2). Small values 
of e correspond to the applicability of the non-relativistic limit. The transverse 
densities of Eq. (|20p results are shown in Fig. Q3 for e = 0.01 and 0.1. As 
expected from non-relativistic intuition, the smaller value of the binding energy 
corresponds to a greater spatial extent. 



2 ZOO OF NUCLEON DISTRIBUTIONS AND DENSITIES 



The vast literature concerning the distribution functions that are used to describe 
nucleon structure includes generalized parton distributions GPDs (|24U25tl2Ul l2"7l 
[MllSSlEQiEIlESlEaiSlESI) and transverse momentum distributions ([36ll371l38l 
39. 40 . 4T|I42|) . The relationship between GPDs, transverse charge densities and 
TMDs is discussed in the present section. 

We begin by discussing generalized parton distributions, GPDs and follow the 
discussion of Burkardt (|14p . Deep- inelastic scattering experiments allow the de- 
termination of parton distribution functions (PDFs), which give the probability 
that quarks carry a given fraction x of the nucleon longitudinal (+)momentum 
in the infinite momentum frame (IMF). PDFs are the forward matrix element of 
a light-like correlation function, 




(20) 



which for the present model is evaluated as 




(21) 




(22) 
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with 

6 q (x,0) = y^I 9 -(-^,0) 7 + g (^,0)e^-, (23) 

where |P, 5) represents the wave function. Here we use the light-cone gauge A + = 
0. The use of the canonical field expansion for the quark field-operators shows 
explicitly that PDFs give the probability that the quarks carry a longitudinal 
momentum fraction x. We do not display the scale dependence of the PDFs to 
simplify the notation. 

Generalized parton distributions (GPDs) (l2ll231l2W7ll2^1l^ 
which describe the scaling limit in real and virtual Compton scattering exper- 
iments, are defined by allowing the momenta and spins of the initial and final 
nucleons to differ: 

{P',S'\6 q (x,0 ± )\P,S) (24) 



. )p+ -u(P , ,s') ^f + H g (x,t,t) + i^ M ^E q (x,t,t) j u(P,s) (25) 

with = \ (P^ 1 + P ,fl ) being the mean momentum of the target, A^ = P'^ — 
the four momentum transfer, and t = A 2 the invariant momentum transfer. 
The skewedness parameter £ = — ^=qr represents the change in the longitudinal 
component of the nucleon momentum. 

GPDs allow for a unified description of a number of hadronic properties. If P' = 
P, the forward limit, they reduce to conventional parton distribution functions 
H q (x, 0,0) = q(x). The integral over x causes the x~ coordinate to vanish, so 
that the operator O q is converted into a local current operator and results in the 
appearance of the usual Dirac form factors of the nucleon 

J2e q JdxH q (x,t,t) = F 1 (t) (26) 
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]Te g fdxE q (x,tt) = F 2 (t). (27) 
q J 

The next step is to define impact parameter parton distributions (|14p . One 



localizes the nucleon in the transverse direction using the superposition Eq. (jlOp 
but with light-cone helicity states ([33]) |P + ,p,A) = \P,S). For a transversely 
localized state |p + ,Rj_ = Oj_,A), Burkardt's impact parameter dependent PDF 
is given by 

q(x,b) = (p+Ri = ± ,\\6 q {x,b)\p + ,-R ± =0 ± ,A). (28) 

The operator 9 (x, b) is obtained by a translation of the operator appearing in 
Eq. ([23]) in the transverse plane 

6 g (x,b) = J^q(-^,b)l + Q(^Me lxp+x ~. (29) 

The function g(x, b) gives the probability density that a quark of momentum 
fraction x is located at a transverse position b. 

Burkardt (|14p used the transversely localized states of the same helicity to 
show that q(x,h) is the two-dimensional Fourier transform of the GPD H q : 

?(*, b) = y e- j qb i?,(x, i = -q 2 ). (30) 

Integration over x using Eq. (|26p and Eq. (|14p shows that 

p(b) = JdxJ2e q q(x,b)= J ^F 1 (t = -q 2 )e- iq ' b (31) 

This means that the transverse density can be obtained either as an integral 
over x of q(x, b) or as an integral over x~ of the three dimensional coordinate 
space density, as in Sect. 1.2. This equality of an integral over a momentum with 
one over a distance is an example of Parseval's theorem. In either case, model- 
independent information regarding the longitudinal coordinate is not available 
because the probe momentum is transverse. That is q + = 0. 
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In the model of Sect. 1.3 (scalar hadron made of two scalar constituents), q(x) 

can be obtained from the integrand of Eq. (|2Q|) as 



x(l — x) 3 



1 — X' 



(32) 



One sees, for large values of x approaching unity, that the transverse extent is very 
narrow because this corresponds to very large values of the relative transverse 
separation B = b/(l — x) for finite values of b. Thus the transverse spatial extent 
depends explicitly on the value of x. In particular, large values of x are associated 
with small values of b. 

Transverse momentum distributions TMDs are another generalization of par- 
ton distribution functions. These, which contain probability distributions re- 
garding both the longitudinal momentum fraction x and transverse momentum 
k carried by the quarks, are given by (|36y 37 [ HU |) 

*J(x = ^,k) = (P,5| J^^^C--^C mrq (C,C)\P,S}, (33) 

where the time separation £ + = 0. Here the quark field operators are the gauge 
invariant operators in gauges for which the gauge potentials vanish at space-time 
infinity (|44|, I45 [) . These operators are necessary for non-zero transverse spatial 
separations, £. The operator T is a Dirac matrix that defines the quark density. 

The initial and final states appearing in Eq. (|33p have the same momentum; 
the major difference between TMDs and GPDs. Another major difference is that 
GPDs involve densities in transverse coordinate space, but TMDs involve densi- 
ties in transverse momentum space. However, one can use Wigner distributions 
(|44p to show that GPDs and TMDs are obtained from different manipulations 
on the same operator. We define the reduced Wigner distribution as 



(C , C, k+ ^) = ^j dv-d 2 V e^'~^q(C - 7f /2, C - r,/2)Tq(C + TA C + f?/2) 
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GPDs are obtained by taking the matrix element of W^{C,~ = 0, C, = 0, k + , k) 

between states of different momentum and then integrating over k as in 

H q (x,£,t) = (P',S\ J ^W^(r,C,k + ,k)\P,S), (35) 

and TMDs are obtained by taking the matrix element of Wq (C~ ; d ^ + ; k) between 
states of the very same momentum: 

*J(x,k) = (P,S\ J ^W?(C,Ck + ,k)\P,S). (36) 

2.0.1 X-SUM RULES, THE CONNECTION TO THE EQUAL TIME FORMALISM, 

and lattice calculations The integral of a GPD over all values of x as in 
Eq. (|26p and Eq. (|27p converts the non-local bilinear appearing in Eq. (|29[) into 
a local operator. More generally, the GPD is a correlation function of quarks at 
the same light cone time x + = 0. The integral over x gives also x~ = 0, so that 
both the time t and spatial z separations vanish. Thus one can compare matrix 
elements involving the integral over x with related matrix elements computed 
using the equal time formalism, t = 0. Moreover, performing the x integral allows 
a connection with lattice QCD calculations. The lattice formulation is built using 
covariant Euclidean space. If one sets it and z to in a lattice calculation one 
obtains a quantity that depends on transverse coordinates suitable for comparison 
with a relevant transverse experimentally measured quantity. This relation has 
been pointed out in connection with transversity observables in Ref. (|46h and 
used in Refs. (|37|H5|H5|I5D|I . 

A similar connection (|5ip between the light-front and equal time formalism 
occurs when using TMDs. A TMD Eq. (|33p gives the probability that a quark 
has a three momentum characterized by (x, k). The relevant matrix element 
involves quarks separated at the same light-cone time £ + = 0. Integration over 
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x sets also Q 
time. 
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to 0, so one obtains a density evaluated using quarks at the same 



3 SINGULAR PIONIC TRANSVERSE CHARGE DENSITY 

Understanding the pion is necessary to learn how QCD describes the interaction 
and existence of elementary particles. As a nearly massless excitation of the 
QCD vacuum with pseudoscalar quantum numbers, the pion plays a central role 
in nuclear and particle physics as the carrier of the longest ranged force between 
nucleons and as a harbinger of spontaneous symmetry breaking. 

Computing the electromagnetic form factor of the pion, F 7T (Q 2 ), at asymptot- 
ically large values of Q 2 from first principles was one of the early challenges to 
using perturbative QCD in exclusive processes ([52 1 HkJIII^H 1551) • Such calculations 
have been extended ((551157115511551) by including effects of higher order in the 
strong coupling constant and higher twist effects. The lowest order results are 
about a factor of three smaller than existing data, and higher order and higher 
twist effects are not small at currently available values of Q 2 ((551) . It is widely 
believed that at large enough values of momentum transfer Q 2 the leading-order 
perturbative formula will be correct, but these large values may be very large 
indeed (|60p . As a result there is considerable experimental interest in determin- 
ing the transition to the region where perturbative QCD can be applied. New 
measurements (]6T[[62p have been performed and more are planned (63). Here we 
review the first analysis (|17p that provided a model-independent pionic transverse 
charge density. 

Recent pion data (]61|j62J provide a very accurate measurement of the longitu- 
dinal part of the electroproduction cross section and the related pion form factor 
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up to a value of Q 2 = 2 Ah GeV 2 . The result is that existing data for the pion 
form factor are well represented by the monopole form 



with R 2 = 0.431 fm 2 . The expression Eq. (|37p allows one to determine the 
transverse density from Eq. (|14p with the result: 



where Kq is modified Bessel function of rank zero. For small values of b this 
function diverges as ~ log(6), hence the transverse density is singular and infinite 
at the origin. This infinity is not to be "cured" by a renormalization procedure 
because the charge density under consideration is the matrix element of a valence 
quark operator between eigenstates of the full Hamiltonian. Furthermore, diver- 
gences of quark distribution functions that occur at small values of Bjorken x do 
not occur here because transverse charge density involves the difference between 
quark and anti-quark densities. Many field theory models, derived even before 
QCD was established (64) obtain a form factor that corresponds to a divergent 
transverse density. Note also that any model, such as vector meson dominance 
or holographic QCD (65,66,67) yielding a monopole form factor has a central 
density with a logarithmic divergence. Thus holographic QCD does not supply a 
representation of the soft component (60) of the pion wave function. 

Intuition regarding a possible singularity in the central charge density may 
be gained from other examples. Suppose that the non-relativistic (NR) limit in 
which the quark masses are heavy is applicable. Then the pion would be a pure 
qq object and the charge density would be the Fourier transform of the form 
factor (Sect. 1.1). Given the form factor of Eq. (|3T|) the NR three-dimensional 



F n (Q 2 ) = 1/(1 + R 2 Q 2 /6) 



(37) 




(38) 
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density would be uniquely given by 

PNrW = 2^ 6 R ' (39) 
where r is the distance relative to the pion center of mass. This density diverges 
at r = 0. If one takes r = \/b 2 + z 2 as demanded by the rotational invariance 
of the non-relativistic wave function, then one finds from Eq. (|15p that the non- 
relativistic transverse charge density takes the form of p(b) of Eq. (|38p . 

The divergence of the central transverse charge density encountered here may 
be the consequence of using a simple parametrization, presently untested by mea- 
surements at larger values of Q 2 . Thus we examine other approaches. Consider 
first perturbative QCD (pQCD), which provides a prediction ([52|) - ([55|) for asymp- 
totically large values of Q 2 that 

lim^MQ 2 ) = l^a s {Q 2 )f 2 JQ\ (40) 

with the pion decay constant f n = 93 MeV, and in leading order: 

with rif the number of quarks of mass smaller than Q and A is a parameter fixed 
by data. The log Q 2 term in the denominator does not lead to a non-singular 
behavior of p(b) for small values of b because p(b) ~ log log r and the pQCD form 
factor corresponds to a singularity at short distance. This singularity would arise 
in any model form factor even one based on sum rules e.g (|60p that approaches 
the pQCD asymptotically. 

Chiral quark models (see the review (|68|) ) present other models (f69l lT0~ |) of 
transverse charge densities that are singular at the center as log b because the 
pion form factor takes the monopole form. It is nevertheless interesting to note 
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that each and every observable quantity, including f n and structure functions, is 

computed to be finite. 

Gaussian models with generalized parton distributions H(x, 0, Q 2 ) (see Eq. ([26])) 
dominated by behavior near x = 1 present a set of examples that also yield a 
form factor with a 1/Q 2 asymptotic behavior, and have a impact parameter 
distribution that is well behaved at each value of x for all b. The key asymp- 
totic features are captured in the simple formula ([TBI ITT]) : H(x, 0, Q 2 ) x ^i = 

(1 - x )n-l e -a(l-xrQ^ with n > 2 SO that q(x, 6)^1 = 2 na(l-x) e ~ b2/iMl ~ Xr) ■ 

This form shows that q(x, b) is well behaved for all values of b and for each value 
of x, but the integral over x contains a logarithmic divergence. 

Relativistic light-front constituent quark models (fT8] PT9 ][T5"]) that describe ex- 
isting form factor data have a non-singular central charge density. These models 
can be most simply derived (|79p by using the impulse approximation (evaluat- 
ing the triangle diagram). One starts by evaluating the integral over the minus 
component of the loop momentum fc^, and then cuts off the remaining integral 
over x = k + /p + , k± using a phenomenological wave function. The wave function 
of Ref. (|75p is chosen to be a power-law form, and the resulting model describes 
the existing form factor data in the space- like region, /„-, and the transition form 
factor f n ry in which a virtual photon transforms a real pion into a real photon. 
The model form factor of (75) and the monopole fit of Eq. (|37p are shown along 
with the measured data in Fig. [2] Both the fit and the model provide a good fit 
to the data, but present very different predictions for larger values of Q 2 where 
measurements remain to be done. The corresponding versions of p(b) are shown 
in Fig. [3] The singularity of Eq. (|38p is manifest as a rapidly rising function as b 
approaches zero, but the relativistic constituent quark model provides a smooth 
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function p(b). The planned experiment f)63[) aims to achieve results up to Q 2 = 6 

GeV 2 . This should be large enough to resolve the differences between the model 

(fT5j) and the monopole fit, or rule both out. However, this will probably not 

be sufficient to reach the regime where pQCD might be valid (|63p . The model 

(fT5j) represents one useful phenomenology, but other interesting models exist. 

Computing the transverse density in those models would be useful. 

4 NUCLEON TRANSVERSE CHARGE DENSITIES 

In previous sections we have emphasized that a proper determination of a nucleon 
charge density requires a relation with the square of a field operator, and that this 
can be achieved through the transverse charge density (|7]li4pi3lli8||g0~|) . p(b), given 
by Eq. (|3ip . We review the first phenomenological analysis of existing data that 
determined p(b) for the neutron and proton (|18|) . The neutron results contradict 
the ancient idea, obtained from both meson-cloud and gluon-exchange models 
(|8iy 82 1[83llM|) regarding the positive value of the non-vanishing charge density at 
the center of the neutron. The meson-cloud idea is that the neutron sometimes 
undergoes a spontaneous quantum fluctuation into heavy proton surrounded by 
a negatively charged pionic cloud, leaving a positive central density (81). The 
quark- model mechanism (82,83,84) involves the repulsive nature of the one gluon 
exchange interaction between two d quarks. 

These well-motivated physical considerations concern statements about the 
three-dimensional Fourier transform of the neutron's electric form factor Gg(Q 2 ). 
However, such a transform is not the charge density. The only model-independent 
charge density is p(b). 

To proceed, recall the definitions of the form factors. With J^(x) as the elec- 
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tromagnetic current operator, in units of the proton charge, the nucleon form 

factors are given by 

(p>,\>\J»(0)\p,\)=u(p',\>) ^ Fl (Q 2 ) + i— q a F 2 (Q 2 )ju(p,X), (42) 

where the momentum transfer q a = p' a — p a is taken as space-like, with q + = 0, 
so that Q 2 = —q 2 = q 2 > 0. The nucleon polarization states are chosen to be 
those of definite light-cone helicities A, A' f)43|) The charge (Dirac) form factor is 
F\, normalized such that -Fi(O) is the nucleon charge, and the magnetic (Pauli) 
form factor is F 2 , normalized such that i^(0) is the anomalous magnetic moment. 
The Sachs form factors ([85]) 

G E {Q 2 ) = Fi{Q 2 ) - ^F 2 (Q 2 ), G M (Q 2 ) = F 1 (Q 2 )+F 2 (Q 2 ), (43) 

were introduced so as to provide an expression for the electron-nucleon cross 
section (in the one photon exchange approximation) that depends on the quan- 
tities G 2 E and G 2 M , but not the product Ge Gm* In the Breit frame, in which 
p = — p', Ge is the nucleon helicity flip matrix element of J°. Furthermore, 
the scattering of neutrons from the electron cloud of atoms measures the deriva- 
tive —dGE(Q 2 )/dQ 2 at Q 2 = 0, widely interpreted as six times the mean-square 
charge radius of the neutron. However, any probability or density interpretation 
of Ge is spoiled by a non-zero value of Q 2 , no matter how small. This is because 
the momentum difference between the initial and final states appears via the use 
of derivatives of momentum-conserving delta functions in the moments computed 
by Ref. (|85p . Any attempt to analytically incorporate relativistic corrections in 
a p 2 /m 2 type of expansion would be doomed, by the presence of the quark mass, 
m q , to be model-dependent. This is explained more thoroughly in Ref. ([86]) . 
We exploit Eq. (|3ip by using measured form factors to determine pip). Recent 
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parametrizations (87, 88, 89, 90"ll9T]) of Ge and Gm are very useful, so we use 
Eq. (|4"3]) to obtain F\ in terms of Ge,Gm- Then p(b) can be expressed as a 
simple integral of known functions: 



A straightforward application of Eq. (|44p to the proton using the parametriza- 
tions of Ref. (I9ip yields the results shown in the upper panel of Fig. [5J The 
curves obtained using the two different parametrizations overlap. Furthermore, 
there seems to be negligible sensitivity to form factors at very high values of Q 2 
that are currently unmeasured. The density is peaked at low values of b, but has 
a long positive tail, suggesting a long-ranged, positively charged pion cloud. 

The neutron results are shown in the lower panel of Fig. |U The curves obtained 
using the two different parametrizations seem to overlap, The surprising result 
is that the central neutron charge density is negative. The values of the integral 
of Eq. (|44|) are somewhat sensitive to the regime 8 < Q 2 < 16 GeV 2 for which 
Ge is as yet unmeasured. About 30% of the value of p(0) arises from this region. 
That p(b = 0) < was confirmed in Refs. (p2l l80ll93ll94"l) . 

The negative central density deserves further explanation. See Fig. 0] in which 
the upper panel shows F\ for the neutron from two parametrizations of Ref. (|9ip . 
In both cases F\ is negative (because of the dominance of the Gm term of Eq. (|44p ) 
for all values of Q 2 . This along with taking 6 = 0, Jo(Qb) = 1 in Eq. (|44p leads 
immediately to the central negative result. The long range structure of the charge 
density is captured by displaying the quantity bp(b) in the lower panel of Fig. |H 
At very large distances from the center, bp(b) < suggesting the existence of 
the long-ranged pion cloud. Thus the neutron transverse charge density has an 




(44) 



with r 



and Jq a cylindrical Bessel function. 
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unusual behavior in which the positive charge density in the middle is sandwiched 

by negative charge densities at the inner and outer reaches of the neutron. A 

simple model in which the neutron fluctuates into a proton and a tt~ parametrized 

to reproduce the negative-definite nature of the neutron's F\ (|95|) reproduces the 

negative transverse central density. In this case, the negative nature arises from 

pions that penetrate to the center. The change from the nominal positive value 

obtained from Ge can be understood as originating in the boost to the infinite 

momentum frame (|86p . 

One can gain information about the individual u and d quark densities by 
invoking charge symmetry (invariance under a rotation by n about the z (charge) 
axis in isospin space) (96,97,98,99) and neglecting effects of ss pairs (100). 

Model independent information about nucleon structure is obtained, with the 
particular surprise that the central density of the neutron is negative. Future 
measurements of neutron electromagnetic form factors could render the present 
results more precise, or potentially modify them considerably. Obtaining a quan- 
titative and intuitive understanding of these results presents a challenge to lattice 
QCD and to model builders. 

4.1 Meaning of central density b = 

The transverse coordinate b = is the transverse center of the nucleon. It is of 
interest to note that, in the infinite momentum frame, this position is also the 
center of the nucleon. In that frame the nucleon has no longitudinal extent. This 
can be seen by considering a light cone wave function ip(x,K). The canonically 
conjugate coordinate to x = k + / P + is P + x~ , (jlOip so that the coordinate-space 
wave function would depend on the product x~P + . Thus when P + is infinite, the 
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coordinate-space wave function will vanish unless x~ = 0. This means that in the 

infinite momentum frame, the longitudinal density appears as a delta function 

and the position b = is the true nucleon center ([6]). 

4.2 Inclusive-exclusive Connection Interpretation of the negative 
neutron charge density 

Generalized parton distributions (GPDs) contain information about the longi- 
tudinal momentum fraction x as well as the transverse position b. Information 
regarding these dependence's is obtained from experiment via GPDs that repro- 
duce both deep inelastic scattering and elastic scattering data. Miller & Arring- 
ton (I102p used this inclusive-exclusive connection to interpret the central neutron 
charge density, with the finding that the center of the neutron is dominated by 
negatively charged d quarks. Their argument is reviewed here. 

The quantities q(x,b) are not measured directly, but have been obtained from 
models that incorporate fits to parton distributions and electromagnetic nucleon 
form factors (|lU3lllU4l[T05l[T06D . Form factor sum rules Eq. ([55|> . Ea. (|27jl at zero 
skewness are exploited to model valence quark GPDs, H% = H q —H q . This yields 
the net contribution to the form factors from quarks and anti-quarks. Possible 
effects of strangeness are neglected in these fits. 

Each parametrization used (103} 11041 1105P incorporates the Drell-Yan-West 
(|1U71I1U8|) relationship between the behavior of the structure function vWi{x) 
function near x = 1, measured in inclusive reactions and the behavior of the elec- 
tromagnetic form factor at large values of Q 2 , measured in the exclusive elastic 
scattering process. In particular, for a system of n + 1 valence quarks, described 
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by a power-law wave function 

]imuW 2 (x) = (l-x) 2n - 1 ^ lim F^Q 2 ) = -L (45) 

The value of n that defines the high-x behavior of the structure function also 
defines the high-Q 2 behavior of the form factor. This relation associates the 
behavior of large values of x with large momentum transfers, Q 2 = q 2 , which 
in turn corresponds to small values of b. There is a further connection between 
large values of x and small values of b which emerges from Eq. (|lip and was seen 
in the example of Sect. 1.3 (and Eq. (|32p ). If a single quark carries momentum 
fraction x very near unity, only one term in the sum of Eq. (jlip survives and the 
restriction of this term to zero, causes the corresponding transverse coordinate to 
vanish. This interesting connection between different components of the position 
and momentum vectors is not associated with the uncertainty principle. 

Given that the d quark dominates the large x quark-distribution function of 
the neutron, see e.g. (|109j) . and that the central charge density of the neutron is 
negative, it is natural to conclude that the central charge density of the neutron 
arises from the predominance of d quarks at the center. This is shown in Fig. 6. 
One sees that for large values of x, b must be small to have a non-zero value of 
p n (x,b), and these non-zero values are negative. 

4.3 Magnetization Density 

We now try use transverse densities to obtain the nucleon magnetization density 
in the infinite momentum frame (110) in terms of a magnetization density. This 
quantity is closely related to the transverse charge density for a polarized nucleon 
obtained by Carlson and Vanderhaeghen (j80|) . Our starting point is the relation 
that fi ■ B is the matrix element of J ■ A in a definite state, \X). 
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Take the rest-frame magnetic field to be a constant vector in the 1 (or b x ) di- 
rection, and the corresponding vector potential as A = Bb y z. Then consider the 
system in a frame in which the plus component of the momentum approaches in- 
finity so that J A —7- J + A~. The anomalous magnetic moment may be extracted 
by taking the matrix element of [i ■ B in the state 

\X) = -J= [\p+, R = 0, +) + \p+, R = 0, -)] , (46) 

where \p + , R = 0, +) represents a transversely localized state of definite P + and 
light-cone helicity. The state \X) (I14p may be interpreted as that of a transversely 
polarized target (lllip . The resulting anomalous magnetic moment \x a is given 
by dUO]) 

» a = (X\ J d 2 bb y ql(0,b)q + (0,b)\X), (47) 

where q(x~,h) is a quark-field operator, and c/+ = 7°7 + g. This matrix element 
is the anomalous magnetic moment because the use of the transversely localized 
state \X) and the factor b y suppresses the Dirac contribution (jTTTj). 

In Ref. (jllOp Ha was evaluated by using Burkardt's (fT4"|) impact parameter 
distribution and then integrating by parts. The result was 

Ma = 2^/ d2bpAl( - b) > (48) 

where 

PM{b) = J = -q 2 )e-^ b , (49) 

and the subscript M denotes the anomalous magnetic moment. 

The expression (|49p has an appealing simplicity. However, the directly-obtained 
Eq. (|47p can also be evaluated immediately. We pursue this here to find 



^=m] dbb »^r> (50) 
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so that the quantity —6, 




jOjvf(b) also has an interpretation as an anoma- 



lous magnetization density. The two integrals appearing in Eq. (|48p and Eq. (|50p 
the same value. Ref. (jllOp rejected the use of /5 a/ (b) as a magnetization density 
because of the appearance of an explicit direction y. However, this direction 
has a general interpretation as the transverse direction orthogonal to that of the 
transverse magnetic field. We evaluate /OA/(b) in terms of F2(Q 2 ) to find 



where <j> is the angle between the direction of b and that of the transverse magnetic 
field, which is also the direction of the nucleon polarization. Thus the physical 
direction of the polarization or magnetic field provides a definite spatial direction. 
Indeed the magnetization density /5a/ (b) is largest in directions perpendicular to 
the direction of the nucleonic polarization (or magnetic field), as shown in Fig. 7. 
The largest values occur for <j> = tt/2, and the magnetization density peaks at 
about 0.5 fm. Furthermore /5a/ (b) vanishes if b = or if = 0. These features are 
in accord with the expectations of classical physics. A current in the z direction 
causes a magnetic dipole density ~ r x J in the x— direction for positions r 
along the y-direction. Therefore we conclude here that the quantity PA/(b) is the 
preferred expression for the magnetization density. 

Hoyer & Kurki have computed the transverse density of the electron (|112p . 
Using their expression for the electron F2 in Eq. (|50p leads to the correct result 
for the electron's anomalous magnetic moment, /t|. Furthermore the presence of 
the sin 2 (j) term is consistent with their physical interpretation of the sign of \i e a as 
caused by currents in the positive (negative) z— direction for positive (negative) 
values of y. Computing the cross product between b and those currents gives 
magnetization in the x-direction for both positive and negative values of y. 




(51) 
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To better understand this new magnetization density, consider the moments of 

both densities. Define the nth moment of pm as 

(b 2n ) M = J d 2 bb 2n p M (b). (52) 

These quantities are related to the nth derivative of F2(q 2 ). Similarly the nth 
moment of pM is given by 

(^)m--/^ 2 ^^- (^) 
Then use of Eq. (j4T)|) and integration by parts shows that 

(b 2n )~ = (n + l)(b 2n ) M . (54) 

The n = moment corresponds to the anomalous magnetic moment which is 
the same for the two densities. The case n = 1 which defines the mean-square 
magnetic radius is more interesting. Ref. (jllOj) showed that (ft 2 ) a/ was slightly 
larger than {b 2 )ch (which is obtained from F\). However {b 2 )^ is twice as large 
as (b 2 )M and therefore is much larger than (b 2 )ch* This clearly shows that the 
proton's magnetization density extends much further than its charge density, a 
conclusion obtained in Ref. (IllOp . This is surprising because it contradicts simple 
naive intuition gained from the rapid fall of the ratio Ge/Gm with increasing 
momentum transfer ([3]), if one assumes that Ge (Gm) is related to the spatial 
extent of the charge (magnetization) density. 

4.4 Non-cylindrically symmetric transverse charge density and 
shape of the nucleon 

Carlson and Vanderhaeghen (I80p computed the transverse charge density p!^ (b) 
in a given state of transverse polarization For the case that the polarization is in 
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the x direction their result is expressed as 



P 

p%{h) =p(b) + sin / 



'OO 



dQ Q 2 
2ir 2M 



Ji(bQ)F 2 (Q 2 ), 



(55) 



where p(b) is obtained from Eq. (131ft . The sign of the second term of Eq. (|55j) 
is obtained using the definition that q is the momentum absorbed by the target 
nucleon. Their results for the proton are shown in Fig. [8j 

The magnetization density p~m(P) presented here and the transverse charge 
density of Ref. (|80|) each depend on the direction of b, thus indicating a non- 
spherical shape of the nucleon or the violation of cylindrical symmetry. These 
violations of cylindrical symmetry should not be confused with those present in 
earlier papers on the shape of the nucleon Refs. (|113yil41l5ip that define a shape 
via the matrix elements of a spin-dependent density operator. The transverse 
charge densities that are the main subject of this article are matrix elements of 
currents that are local in transverse coordinate space and are related to GPDs. 
In contrast, those shapes of (|113pil4"ll5ip . computed from momentum-space wave 
functions and probabilities are matrix elements of quark densities for given values 
of transverse momenta and are related to TMDs. It should also be noted that the 
first example of (|113p is presented in coordinate space. The infinite momentum 
frame version of the spin-dependent density operator is q\_(0, b)i(l+n-7)g + (0, b), 
where n is an arbitrary transverse direction f)49|) . Evaluating the matrix element 
of this quantity produces the spin-dependent density, which can be thought of as 
the x~ integrated version of the coordinate space results of (113). The proton is 
non-spherical if the function A'^ 10 of (|50p is non-vanishing, as recent lattice calcu- 
lations show (|115p . One can also examine several different generalized densities 
q+(0, b)Tq + (0, b), with operators T denoted in (j4"9j) . 
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5 TRANSVERSE TRANSITION CHARGE DENSITIES 

There is much experimental interest in measuring transition form factors that 
represent probability amplitudes for a nucleon undergoing electron scattering to 
make a transition to a given baryon resonance. Substantial experimental efforts 
are underway at electron-scattering laboratories such as Jefferson Laboratory, 
ELSA in Bonn, and MAMI in Mainz. Moreover, there is considerable interest in 
using lattice techniques to make QCD calculations of these transition amplitudes. 

Carlson & Vanderhaeghen (|80p used empirical information regarding the N — > 
A transition form factors to map out the transition charge density. This, in 
a transversely polarized N and A, contains monopole, dipole and quadrupole 
patterns. The latter corresponds to a deformation of the nucleon and A transverse 
charge density. Substantial deformations are observed. Lattice QCD calculations 
have been applied to these very same transition densities (|47p48|) . The A + charge 
density is found to be elongated along the axis of the spin and the quadrupole 
moment is larger than the value characterizing a point particle. This means the 
the A + is prolately deformed. 

Tiator & Vanderhaeghen (|116p used recent experimental data to analyze the 
electromagnetic transition from the nucleon to the Pn(1440) resonance, which 
is often believed to be the first radial excitation of the proton. They used the 
empirical transition form factors to find that the transition from the proton to 
the Pn(1440) is dominated by up quarks in a central region of width of about 
0.5 fm and by down quarks in an outer band which extends up to about 1 fm. 
Tiator et al. (jTTTJ) are extending the study of transition transverse densities to 
the Sn and D13 resonances. 
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6 NUCLEAR TRANSVERSE CHARGE DENSITIES 

The very same infinite momentum frame formalism can be applied to the analysis 
of nuclear form factors. The non-relativistic approximation of Eq. (J6|) is generally 
applicable for existing measurements of nuclear form factors ©. However, meson 
exchange currents are known to be important for Q 2 > 2 GeV 2 (|118p . In that 
case, the charged constituents have different masses (see Eq. 0), the pion moves 
relativistically, and the charge density is not a three-dimensional Fourier trans- 
form of the electromagnetic form factor. Moreover, Jefferson Laboratory plans to 
considerably extend the available range of momentum transfer for light nuclear 
form factors. Hence transverse charge densities which supply model independent 
information are of considerable relevance for nuclear physics. 

Carlson & Vanderhaeghen (|119p made the first computation of a nuclear trans- 
verse density in their analysis of the deuteron empirical transverse charge densities 
p2U|) . The charge densities are characterized by monopole, dipole and quadrupole 
patterns because the deuteron has a spin of unity. Ref. (|119p observes a dip in the 
center-of-charge density for helicity-zero deuterons. This central depression is in 
accord with standard nuclear force model calculations (|121 1) and is a simple con- 
sequence of the deuteron D-state. We note that the deuteron is a non-relativistic 
system in the sense defined by Ref. (|6j). Thus these densities can, with a high 
degree of accuracy, be thought of as an integral over z of the corresponding non- 
relativistic density; see Eq. f)15|) . The dip is an important finding, because it is 
model independent. Ref. (|119p also finds that transversely polarized deuterons 
show dipole and quadrupole structure in the charge densities. Their electric 
dipole and quadrupole moments only depend on the spin-1 particle's anomalous 
magnetic dipole moment and its anomalous electric quadrupole moment, arising 
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from its internal structure. 

Nuclei are non-relativistic systems ([6]) so that for presently available data the 
principle distinction between transverse densities and the usual three-dimensional 
Fourier transform is that the former is an integral Eq. (|15p over the latter. How- 
ever, existing nuclear data extend only to about Q 2 = 1 GeV 2 (|122yi23pi24p . In 
this case the motion of the nucleons within the nucleus is just barely relativistic. 
If, as expected, data are taken at higher values of Q 2 relativistic effects can be 
expected to be important. One example is the nuclei with ^4 = 3 and mass M3 
and Sachs electric and magnetic form factors Fc,Fm (122). If Q 2 /AM 2 is of 
order unity, the distinction between F\ = Fq — Q 2 /4:M 2 Fm and Fq will become 
important and the role of transverse densities may become very important. 

6.1 Other Applications 

Transverse momentum densities, giving the momentum density of hadrons, have 
also been studied (|125[|126|) . The spatial distribution of the momentum com- 
ponent P + is found to be related to Fourier transforms of gravitational form 
factors (|125|) which in turn are related to experimentally observed data. This is 
because gravitational form factors are matrix elements of the energy-momentum 
tensor obtainable as second Mellin moments of GPDs (|127^ 25 ^128|) . The proton 
momentum density in the transverse plane plane was found to be more compact 
than its charge density (125). 

Any charge density will deform when subjected to an external electric field; 
related measurable quantities are called polarizabilities. Gorchtein et al. (|129p 
extended the transverse charge density formalism to extract light-front quark 
charge densities related to polarizabilities and showed that the resulting induced 
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polarizations can be extracted from proton generalized polarizabilities (1301 1131]) . 

The available data for the generalized electric polarizabilities of the proton yield 

a transverse density that shows an usual oscillatory pattern. 

7 SUMMARY AND FUTURE DIRECTIONS 

Transverse charge densities are a new tool for analyzing electromagnetic form 
factors of systems composed of constituents that move relativistically. These 
quantities are hadronic charge densities as seen in a reference frame moving 
with infinite momentum. One of the main advantages of the infinite momen- 
tum framework is that boosts in the transverse direction form a kinematic sub- 
group of the Poincare group, so that transformations to frames moving in a 
direction transverse to that of the infinite momentum are carried out using the 
transverse position operator. This is just like the usual non-relativistic Galilean 
transformation. Thus the two-dimensional Fourier transformation of electromag- 
netic form factors provides a rigorous way to study charge and magnetization 
densities. Transverse charge densities involve matrix elements of local operators 
q\_(t = 0, z = 0, h)Tq + (t = 0, z = 0, b), computable using any of the light-front, 
equal-time or lattice techniques. 

The use of transverse momentum densities has led to some very interesting 
findings. Examples include the negative nature of the central neutron charge 
density (Tl8|) as caused by the dominance of d quarks at the center fll02|) , the 
seemingly singular nature of the density at the center of the pion (fT7|) , the spatial 
extent of the proton's magnetization density is greater than that of its charge 
density (jllOp . and the deformation of the nucleon fjl 13^11 15]) and the A baryon 
(|8U|) is substantial. The positive sign of the electron anomalous magnetic moment 
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is explained (I112p . 

Transverse densities have been used to analyze other baryon transitions (116), 
nuclear charge densities (|119p momentum densities (|128p and generalized polar- 
izabilities (fT29|) . 

Measurements of form factors at larger values of Q 2 than are presently avail- 
able is needed to test the singularity of the pionic transverse density and the 
inner core of the neutron transverse density. It also seems likely that the use 
of transverse charge densities will become the chosen tool to analyze nuclear 
charge distributions once very high-momentum transfer data become available. 
Transverse densities are relevant whenever matrix elements of an operator q\(t = 
0,z = 0,h)Tq + (t = 0, z = 0, b) appears. Thus one may expect to see a host of 
interesting, informative and important future applications of transverse densities. 
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Figure 1: bp{b) for the model of Eq. ([2D]). Solid e = 0.01,dashed e = 0.1. 
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Figure 2: Pion form factor data as plotted in (|62p. The data labeled Jlab are 
from (62). The data Brauel et al. 1171 J) and that of Ackermann et al. (j72j) have 
using the method of (62). The Amendola data et al. are from (73) The data 
point labeled PionCT is from (|74p . The (red) dashed curve uses the monopole 
fit Eq. (I37p and the (black) solid line the constituent quark model of (|75h . bp(b) 
for the model of Eq. (120j) . Reprinted from Ref. (|17h with permission of the APS. 
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Figure 3: p(b) corresponding to the two models shown in Fig. [2j The (red) 
dashed curve represents the transverse density obtained from the monopole fit 
and the (blue) solid-line is obtained using the relativistic constituent quark model 
of ([75]) . Reprinted from Ref. (fl~7|) with permission of the APS. 
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Figure 4: Nucleon p(b) Upper panel: proton transverse charge density. Lower 
panel: neutron transverse charge density. These densities are obtained using the 
parametrization of (|9ip . 
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Figure 5: Neutron F\ and bp(b). Upper panel: F\{Q 2 ). Lower panel: bp(b) . The 

solid curves are obtained using Fit 1 and and the dashed curves with Fit 2 of 

p„(x,b) 




Figure 6: Neutron impact parameter charge distributions, p n (x,b) as a function 
of b. 
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Figure 7: Upper panel: /?A/(b = (b, <f> = vr/2)) as a function of b. Lower panel: 
Density plot of /5m (b) • The horizontal axis is the direction of the applied mag- 
netic field. The largest (smallest) values of pu are denoted by the brightest 
(darkest) areas. This figure is obtained using a dipole parametrization for F<i of 
the proton. 
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Figure 8: Quark transverse charge densities in the proton, after (|80p . figure 
provided courtesy of M. Vanderhaeghen. The upper panel shows the density in 
the transverse plane for a proton polarized along the x-axis. The light (dark) 
regions correspond with largest (smallest) values of the density. The lower panel 
compares the density along the y-axis for an unpolarized proton (dashed curve), 
and for a proton polarized along the x-axis (solid curve). The proton form factors 
are from Arrington et al. (|9Up . The momentum transfer qis the momentum added 
to the target proton. 



